Accelerating f(T) gravity models constrained by recent cosmological data 
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Generalised Teleparallel gravity, also referred to as f(T) gravity, has been recently proposed as 
an extended theory of gravitation able to give rise to an accelerated expansion in a matter only 
universe. The cosmic speed up is driven by an effective torsion fluid whose equation of state depend 
on the f(T) function entering the modified gravity Lagrangian. We focus on two particular choices 
for f(T) which share the nice property to emulate a phantom divide crossing as suggested by 
some recent data. We check their viability contrasting the predicted background dynamics to the 
Hubble diagram as traced by both Type la Supernovae (SNela) and Gamma Ray Bursts (GRBs), 
the measurement of the rate expansion H(z), the Baryon Acoustic Oscillations (BAOs) at different 
redshifts, and the Cosmic Microwave Background Radiation (CMBR) distance priors. Both f(T) 
models turn out to be in very good agreement with this large dataset so that we also investigate 
whether it is possible to discriminate among them relying on the different growth factors. 
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I. INTRODUCTION 



The discovery of the acceleration of the universe 
through the SNela Hubble diagram [l[ , recently awarded 
by the Nobel prize, has been latter confirmed by wide 
range of data, from more recent SNela data to BAOs and 
CMBR anisotropies @, Q . On the other hand, such over- 
whelming abundance of observational evidences in favour 
of the cosmic speed up does not fit in the framework of 
General Relativity (GR) making clear that out theoreti- 
cal background is seriously flawed. The naive interpreta- 
tion that an unexpected new ingredient in the form of a 
negative pressured component is driving the acceleration 
poses difficult questions on its nature and nurture and 
introduces further problems hard to be solved. It has 
therefore gained more and more attraction the hypothe- 
sis that the cosmic speed up is rather the first signal of 
a breakdown of our understanding of the laws of gravity 
on cosmological scales. Motivated by this consideration, 
much attention has been dedicated to f(R) theories Q 
where the scalar curvature R is replaced by a suitably 
chosen function f(R) in the gravity Lagrangian. 

Looking for a correction to GR, it is instructive to re- 
member that an equivalent formulation is represented 
by teleparallelism. In this theory, torsion, instead of 
curvature, is responsible of the gravitational interaction 
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0, H| and the Weitzenbock connection replaces the Levi - 
Civita one on the underlying Riemann - Cartan space- 
time. In this scenario, gravitational interaction is not 
replaced by geometry and the torsion acts as a force, al- 
lowing the interpretation of gravity as a gauge theory of 
the translation group Q- Despite conceptual differences, 
teleparallel gravity and GR yield equivalent dynamics so 
that the interpretation of the gravitational interaction in 
terms of a curved or torsioned spacetimc is only a matter 
of convenience, at least at the classical level. 

In as much the same way as for f(R) theories, one can 
obtain a generalised teleparallel gravity replacing T with 
a generic function f(T) thus opening the way to a rich 
phenomenology. A particular important consequence is 
the breakdown of the equivalence with the classical GR 
with the two theories now predicting a radically different 
dynamics [||. Modified teleparallel gravity preserves the 
advantage of giving equations that are still second or- 
der in field derivatives opposite to the fourth order equa- 
tions deduced in f(R) gravity thus avoiding unpleasant 
pathologies. On the other hand, it suffers from the lack 
of Local Lorentz Invariance (LLI) so that all the 16 com- 
ponents of the vierbien are independent and one cannot 
simply fix 6 of them by a gauge choice [9|. 

A critical role in generalised teleparallel theories is 
played by the choice of the functional expression for 
f(T). The lack of firmly established theoretical con- 
straints leaves open the way to a wide range of possi- 
bilities which can only be validated a posteriori, i.e. by 
contrasting their predictions with the observational data. 
This is the aim of the present work where we focus our 
attention on two particular classes able to give rise to a 
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phantom - like behaviour of the effective torsion fluid. We 
then test these two models against SNela + GRB Hubble 
diagram, H(z) measurements from cosmic chronometers, 
BAOs data and the CMBR distance priors. Although 
wide, the present dataset only traces the background ex- 
pansion so that we will also investigate whether further 
insight into the properties of these models can be ob- 
tained by the analysis of the growth factor being this 
latter a quick way to look at how perturbations evolve in 
the proposed modified tclcparallel scenarios. 

The layout of the paper is as follows. In Scct.HH we 
briefly review the cosmology of f(T) gravity and present 
the two models we are going to analyse. Sect. HIT] is de- 
voted to the likelihood analysis and the results, while 
Sect. lIVI enlarges the study to the growth factor. Finally, 
conclusions and perspectives are drawn in Scct.fVl 



II. f(T) GRAVITY 

Teleparallelism is a dynamical theory for the vierbein 
{e a }, whose components in a given coordinate basis 
help in introducing the metric as a subsidiary field : 
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Vabe a Jx)e b l/ (x) , 



where n a i, = diag(l, — 1, — 1, — 1). The dynamics is then 
described by the action 



S = TTTT^ / [T + f(T)]ed 4 x + S M 
167rG J 



(1) 



where e = det e p = ^/—det(gp U ) and Sm is the action for 
the matter fields. In Eq.([T]), a key role is played by f(T) 
which is a diffcrcntiablc function of the torsion tensor T 
defined as 



rp Qp rjipu 

^ pv p 



One can naively expect that, in order to get the cos- 
mological counterpart of the field equations, one has sim- 
ply to insert the Robertson - Walker (RW) metric into 
Eqs.®. Actually, some care is needed since, due to the 
lack of LLI, two pairs of vierbein that lead to the same 
metric tensor are not equivalent from the point of view 
of the theory. Nevertheless, in case of spatially flat RW 
metric, a convenient choice is represented by the diagonal 
vierbein, i.e. fill 



dt 



a(t)dx l 



where a(t) is the scale factor as function of cosmic time 
t. With such a choice, the dynamical equations become : 



H 2 = ^ P -\f{T)-2H*f T {T) 



(H2 y = 16nG P + 6H 2 + f(T) + 12H 2 f T (T) 



(3) 



2AH 2 f TT {T) - 2 - 2/ T (T) 



while the torsion scalar reduces to T = —6H 2 with H = 
a/ a the usual Hubble parameter and (p,p) the energy 
density and pressure of the matter component. Note that 
hereafter we will denote with a prime and with a dot 
differentiation with respect to In a and t, respectively. 

The modified Fricdmann equations ([3]) can be rewrit- 
ten in the usual form by introducing an effective dark 
torsion fluid with energy density reading 



Pt 



2T/ t (T) - f(T) 
16ttG 



(4) 



with 



Since matter still minimally couples to gravity only, its 
conservation equation will be unaffected so that we still 
have pm oc a -3 and p r oc a~ 4 for the scaling laws of 
dust matter and radiation. Imposing the Bianchi identi- 
ties and the diffcomorphism invariance of the theory, the 
conservation equation for the effective torsion fluid reads 



I, 1 1 Pt + 3if(l + w T )p T = 

QP — _ (TP —TP 4- TP \ J f)PT a t) p T a 

LJ [IV A V p,U pv ~ up) ' ^ p CTf rj V 1 <jp 1 

having set 



Varying the action with respect to the vierbein compo- 
nents e®(x), one gets the field equations 

e~%(e e» a S p mi + h) + e x a S p ^ x (l + f T ) 
+ef a S p ^d^T)f TT + h:(T + /) = 47rGe^r;, (2) 

where T" is the matter energy - momentum tensor and 
subscripts T denote diffcrentation with respect to T. 



= f/T -Jt + 2Tf TT + (n r /3)(/ T + 2Tf TT ) 
WT (1+/t + 2T/ tt )(//T-2/ t ) U 

for the equation of state (eos) of the dark torsion fluid. 
Note the coupling to the radiation energy density through 
the term J7 r (a) = 8nGp r {a)/3H 2 (a). For f(T) = 0, one 
has pr = and tclcparallel gravity goes back to the stan- 
dard GR, while the choice f(T) = const gives wt = — 1 
and the ACDM model is recovered. 

Moving away from the naive choice f(T) — const gives 
rise to rich phenomenology in the dynamical evolution of 
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FIG. 1: Present day value of the torsion eos for the tanh model 
setting Q M = 0.20 (dot dashed blue), 0.25 (solid black), 0.30 
(dashed red) and the fiducial fi r used in the text. 



generalised teleparallel theories. There are actually al- 
most no theoretical hints on the functional form of /(T) 
with, on the contrary, many possible expressions leading 
to an accelerated expansion. It is particularly interesting 
to look at models which are able to give rise to an effec- 
tive cos (defined later) crossing the phantom divide, i.e., 
w e ff(z) > — 1 for z > Zp/j with z p h the phantom divide 
redshift. Two recently proposed model of this kind can 
be obtained setting [12| 



( y 



f(T) = 



a(-T)> 



To 

- exp ( -p— 



(6) 



where the subscript denotes present day quantities. 

In the first case, which we will refer to as the tanh 
model, one should set n > 3/2 in order to have a dark 
torsion fluid with a positive energy density, while the 
same condition leads to n > 1/2 for the second choice, 
referred to hereafter as the exp model. The constant a 
can be expressed as a function of the other parameters 
by inserting Eq. (|6]) into the Friedmann equations ([3]) and 
imposing E 2 (z = 0) = H 2 (z = 0)/H$ = 1. This gives: 



(6# 2 r n (i 



n 



M 



n r ) 



2sccli (£ 2 = !) + (!- 2n) tanh (E 2 = 1) 



(1 — — 



1 - 2n - (1 - 2n + 2p)cP 



(7) 



FIG. 2: Contour plots in the (n,p) plane imposing Wt(z = 
0) = —1.0 for the exp model with 51a/ = 0.20 (dot dashed 
blue), 0.25 (solid black), 0.30 (dashed red) and the fiducial 
Q r used in the text. Note that the lower set of contours are 
almost perfectly superimposed. 



for the tanh and exp model. Note that we have dropped 
the subscripts for the dust and radiation present day 
density parameters (il^j,fl r ) to simplify the notation. 

It is instructive to look at the effective torsion fluid eos 
Wt{z) for the two f(T) models we are considering. We 
do not give here the lengthy expressions which can be 
obtained by inserting Eqs.© into the definition ([5]), but 
show how its present day value wt(z = 0) depend on the 
f(T) parameters in Figs.Q] and [2] for the tanh and exp 
model, respectively. Since the ACDM model is known to 
very well fit the data, it is reasonable to expect that, in 
order to fit the same data, the parameter space of both 
models will collapse in a region giving wt(z = 0) = — 1. 
Fig-ffl shows that, for the tanh model, this will favour 
values of n ~ 1.6 and that the larger is n, the smaller 
fijvf must be in order to have a today effective eos close 
to the ACDM one. This same requirement works in a 
different way for the exp model. Fig.[5J indeed, shows 
that, for given (Omj n ) values, there can be more models 
with wt{z = 0) = — 1 but different p. Put in other words, 
the (n,p) parameters are intrinsically degenerate so that 
the condition wt(z = 0) = — 1 is unable to discriminate 
among different parameter sets. Actually, when fitting 
the data, we do not only look at the present day w e ff 
value, but also to its variation with z. As we will see, this 
helps breaking (at least partially) the degeneracy among 
(ji,D,m) for the tanh model and (n,p) for the exp one. 
Figs.[T]and[2]arc nevertheless useful to quickly grasp what 
is the region of the parameter space of greater interest for 
both the tanh and exp f(T) models. 
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III. f(T) MODELS VS DATA 



The above considerations on the effective eos suggest 
that the proposed generalised tclcparallel scenarios can 
represent viable alternatives to the ACDM model. To 
check whether this is indeed the case, we contrast them 
with a wide set of cosmological data as described below. 



A. Likelihood analysis 

In order to answer the question whether f(T) gravity 
can reproduce the observed Universe, we will explore the 
model parameter spaces by investigating the following 
likelihood function 



£(p) = £ (ti (p) x £ff(p) x C B ao{p) x £ C mb(p) (8) 
where the set of model parameters is 



P = 



(ClM,h,n) 
(fiM, h, n,p) 



for the tanh and exp models, respectively, h being the 
Hubble constant Hq in units of 100 km/s/Mpc, and we 
have set the present day radiation density parameter as 



y h~ 2 (l + 0.2271 AT e// ) 



with (w 7 , N ef f) = (2.469 x 10" 5 , 3.04) in agreement with 
WMAP7 3 constraints. All the terms in the full likeli- 
hood can be written as : 



A(p) 



exp [-X?(p)/2] 
(27r)^/ 2 r 4 1/2 



where (J\fi,T{,Xi) depend on the dataset used. 

The first term refers to the Hubble diagram, i.e. the 
distance modulus \x as function of the redshift z. This is 
related to the underlying cosmological model as 

H(z) = 25 + 5 log d L (z) = 25 + 5 log [(1 + z)r(z)] (9) 
with (1l(z) = (1 + z)r(z) the luminosity distance and 



r ( z ) =— r 

Ho Jo 



dz' 



E(z',p) 



(10) 



the comoving distance. As a tracer, we use both the 
Union2 SNela dataset [lU, comprising AfsNeia = 557 
objects wiht 0.015 < z < 1.40, and N G rb = 64 GRBs 
probing the redshift range (1.48,5.60) with the distance 
modulus estimated in 14 1 starting from the data in [l5| . 
For both dataset s, it is 



Mi 

n 



AG 

E 



where aj is the measurement error for the j-th object, 
while ai n t takes care of the intrinsic scatter of the tracer 
around the relations used to estimate its fi value. For 
SNela, we set ai n t — 0, while, for GRBs, this is a nuisance 
parameter we marginalize over. 

While the Hubble diagram probes the integrated ex- 
pansion rate, the second dataset we use explicitly refers 
to the Hubble parameter H(z) as determined through 
the differential age method pjf. By using red envelope 
galaxies as cosmic chronometers [13], Stern ct al. have 
estimated H{z) over the redshift range 0.10 < z < 1.75 
[l8l ]. To this sample, we add the Ho determination from 
local distance ladders derived by the SHOES collabora- 
tion [l{| and define a total % 2 as for the SNela one. 

While SNela and GRBs probe the distance - redshift 
relation as standardizeable candles, Baryon Acoustic Os- 
cillations (BAOs) work as standard rulers. We therefore 
add the term Lb AO to the full likelihood following the 
method detailed in (20|. To this end, three different sets 
of data are used. First, one use the 6dFGRS [2l| and 
the SDSS [lU surveys to determine the scaled volume 
distance parameter 



r s {z d ) 
dv(z) 



r s (z<i) x 



czr 2 (z) 



H E(z) 



(11) 



with the sound horizon to distance z given by 



.(*) 



E- 1 {z)dz' 



v /l + (3w b )/(4w r )(l + z')- 1 



(12) 



and Zd is the drag redshift. The d, value at z = 0.106 and 
its error is taken from [2l|, while [22| gives d z for z = 0.20 
and z = 0.35 with the corresponding covariance matrix. 
Further BAOs constraints come from the WigglcZ survey 
[20| which recommend to use as observable quantity the 
acoustic parameter [23| 



A{z) 



y/ShfH$dv(z) 



cz 



(13) 



with the volume distance dy{z) given in Eq. ([TTj) . We use 
the observed values and their covariance matrix for A(z) 
determinations at z = (0.44,0.60,0.70) reported in [201 ] . 
We finally define the BAO likelihood as 
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Id xbf (x) x 68% CL 95% CL 



Id xbf (x) x 68% CL 95% CL 



n M 0.286 0.286 0.287 (0.274, 0.299) (0.264, 0.311) 
h 0.719 0.722 0.722 (0.712, 0.734) (0.702, 0.745) 
n 1.616 1.610 1.615 (1.581, 1.636) (1.547, 1.667) 



TABLE I: Constraints on the parameters for the tanh model. 
Columns are as follows : 1.) id parameter, 2.) best fit, 3.) 
mean, 4.) median, 5.), 6.) 68% and 95% confidence levels. 



Qm 0.284 0.286 0.287 (0.276, 0.297) (0.265, 0.308) 

h 0.724 0.731 0.731 (0.723, 0.740) (0.713, 0.749) 

n 1.152 0.757 0.736 (0.577, 0.939) (0.514, 1.103) 

p 0.814 -0.110 -0.100 (-0.263, 0.046) (-0.395, 0.131) 

TABLE II: Same as Table U but for the exp model. 



with 
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^SDSS 
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where D is a Mi dimensional vector with the difference 
between observed and predicted values and Ci the cor- 
responding covariance matrix. 

The last term in the likelihood finally refers to the 
WMAP7 distance priors which have been recommended 
as a quick and efficient way to include the CMBR con- 
straints without computing the full anisotropy spectrum. 
Following 0], we then define the CMBR likelihood sim- 
ilar to the SDSS and WiggleZ ones above, but now the 
observable quantities are the redshift to the last scat- 
tering surface, computed using the approximated formu- 
lae in [24j . the acoustic scale I a = nr(z+)/r s (z+) and the 
shift parameter 1Z [25| 



either as an actual fluid or an effective one induced by 
nonlinear terms in the gravity Lagrangian. Although the 
additional f(T) term does not vanish at high z for all 
the possible combinations of the model parameters, it is 
nevertheless easy to check that the effective torsion fluid 
energy density is by orders of magnitude smaller than the 
matter one at z* so that we can still rely on the WMAP7 
determination of the distance priors. 

In order to efficiently explore the parameter space, we 
use a Markov Chain Monte Carlo (MCMC) code running 
multiple chains and checking the convergence through the 
Gclman - Rubin criterium [26| . The best fit parameters 
will be the ones maximizing the full likelihood, but the 
most reliable constraints on each single parameter pi are 
obtained by marginalizing over all the parameters but the 
i - th one. We will give the mean and median as reference 
value and use the 68 and 95% confidence ranges as our 
1 and 2a errors (although this is formally correct only if 
the marginalized likelihood is a Gaussian function). 



B. Results 

The best fit parameters and marginalized constraints 
for the tanh and exp model are given in Tables U and UH 
while Fig. |3] shows the remarkable agreement among the 
best fit models predictions and the SNela + GRB Hubble 
diagram and H(z) data. The overall good quality of the 
fit may be further appreciated by comparing the model 
predictions for the BAO and CMB quantities with the 
observed values. For the tanh model, we get 



n = 



c/H 



(14) 



Note that, in order to use the distance priors, one should 
implicitly assume that the early universe is matter dom- 
inated with no significant contribution from dark energy, 



f d b f W6 = 0.3418 vs cffi^nR = 0.336 ± 0.015 



0.106 



d b f 200 = 0.1864 vs dg h | 00 = 0.1905 ± 0.0061 , 



. d O350 = °- 1173 vs d o%o = 0.1097 ± 0.0036 
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FIG. 3: Best fit curves superimposed to the SNela + GRB Hubble diagram (left) and H(z) data (right) for the tanh (red solid) 
and exp (blue dashed) models. Note that the fi(z) curves are almost perfectly superimposed so no difference is seen in the plot. 



' A bf (0.U) = 0.467 vs A obs (0.U) = 0.474 ± 0.034 

.4 b/ (0.60) = 0.442 vs .4 obs (0.60) = 0.442 ± 0.020 

t ^4 6/ (0.73) = 0.422 vs A obs (0.73) = 0.424 ±0.021 

l h [ = 302.66 vs if s = 302.09 ± 0.76 

K bf = 1.733 vs Tl obs = 1.725 ± 0.018 , 



z b J = 1092.04 vs z obs = 1091.30 ± 0.91 



so that the best fit tanh model predictions are well within 
lcr from the observed values. For the exp model, we get 



jbf 



0.3428 vs dg 6 f 06 = 0.336 ± 0.015 



d b J om = 0.1865 vs d%% = 0.1905 ± 0.0061 



> d a350 = °- 1121 vs d o%o = °- 1097 ±0.0036 

' _4 b /(0.44) = 0.465 vs A obs (0A4) = 0.474 ± 0.034 

_4 b/ (0.60) = 0.439 vs .4 obs (0.60) = 0.442 ± 0.020 

t A bf {0.7S) = 0.418 vs „4 obs (0.73) = 0.424 ±0.021 

l h [ = 302.54 vs if s = 302.09 ± 0.76 
K bf = 1.735 vs K obs = 1.725 ± 0.018 , 



1092.12 vs z 



obs 



1091.30 ±0.91 



still in good agreement with the BAO and CMBR data. 

The left panel in Fig. [3] is quite instructive showing us 
how the two best fit models can not be discriminated at 
all based on the Hubble diagram data. This can be easily 
understood looking at the right panel where one can see 
that the rate expansion H (z) is quite similar among the 
two models with the largest deviations occurring in the 
rcdshift range (0.5, 1.5). Being of order only few %, it is 
easy to understand that the integration needed to com- 
pute the distance modulus easily washes out the differ- 
ences between the two models thus leading to the almost 
perfectly superimposed curves shown in Fig.[3J 

Another way to explain why the two models make 
so similar predictions is to look at the reconstructed 
effective torsion fluid eos. Fig. [4] shows wt(z) as con- 
strained by the likelihood analysis. In each panel, the 
red solid curve refers to the best fit model, while the blue 
dot -dashed (dashed) one gives, for each z, the median 
value (the 68% confidence range) as estimated evaluating 
wt(z,p) for the model parameters in the MCMC chain. 
For both best fit models (i.e., look at the red curves), the 
eos stays close to the cosmological constant value w = — 1 
up to z ~ 1 so that the dynamics is approximately the 
same as the ACDM one. On the other hand, wt(z) has 
a radical different behaviour for the two models in the 
higher rcdshift regime. Indeed, while the tanh effective 
torsion fluid has an increasing eos, the exp model eos 
stays roughly constant to a value slightly smaller than 
the present day one thus better mimicking a ACDM sce- 
nario. Nevertheless, both best fit models make similar 
predictions for the CMBR distance priors as a result of 
the very small contribution of the torsion fluid to the 
energy budget at the last scattering surface. 

As a final remark, it is worth spending some more 
words on an unusual feature of the exp model. As one 
can see from TableHTl the best fit (n,p) parameters are 
strongly different from their median values with p being 
outside the 95% confidence range. As a consequence, the 
reconstructed torsion fluid cos for the best fit model is 
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FIG. 4: Constrained effective torsion fluid eos wt{z) for the tanh (left) and exp (right) f(T) models. 



outside the 68% confidence range delimited by the blue 
dashed lines in the right panel of Fig. [4] Such an unex- 
pected result is actually a consequence of the strong de- 
generacy among the model parameters. As already high- 
lighted when commenting Fig.[5J for a given n value, there 
are more than one p value giving the same wt{z). This 
intrinsic degeneracy can only partially be broken by fit- 
ting data at different redshifts. However, since most part 
of the data probes the range z < 1, all the couples (n,p) 
giving wt(z) ~ — 1 over this range arc allowed by the fit. 
This is indeed the case for the exp models with (n, p) set 
to their best fit values (leading to wt{z = 0) = —1.042) 
or to the median ones (giving wt(z = 0) = —1.038). For 
larger z, the median and best fit models eos are radi- 
cally different, but the predicted distance priors values 
are nevertheless almost the same since the torsion fluid 
becomes subdominant so that the details of its eos turn 
out to be meaningless. In the following, we will however 
set the {n,p) parameters to their median values since 
such a choice more reliably takes into account how the 
data constrain the exp model parameters space. 



IV. THE GROWTH FACTOR 

The above analysis, while confirming that both /(T) 
models nicely work at reproducing the available data, has 
also highlighted that probing only the background expan- 
sion does not allow us to conclusively discriminate be- 
tween them and the ACDM scenario. In order to deepen 
the investigation and look for a way to put stronger con- 
straints, one has to resort to the analysis of the evolution 
of the perturbations. As a first step, we consider here 
the growth factor g = din 5/ din a with S = Spm/pm the 
matter density contrast. The 6 evolution can easily be 
determined solving [27} 

S m + 2H6 M - ^G e}f p M 5M = , (15) 
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FIG. 5: Growth factor g(z) for the tanh (red solid) and exp 
(blue dashed) f(T) models superimposed to the [29l. |30| data. 



which is formally identical to the ACDM one provided 
the Newton constant G is replaced by the effective one 
G e ff = G/(l + /V). It is worth stressing that Eq. ([T5|) ac- 
tually holds only in the sub -horizon limit, i.e., under the 
approximation k >> /%,'H 2 with T-L the comoving 

expansion rate. A much more complicated expression |28| 
holds in the general case and must be used in evaluat- 
ing large scale probes as the matter power spectrum and, 
for instance, weak lensing obscrvablcs. However, for this 
preliminary investigation, we will use the simpler Eq. (|15[) 
which is well suited for galactic and galaxy cluster scales. 

Fig. [5] shows the growth factor for the best fit tanh and 
the median exp models together with a compilation of 
g(z) measurement taken from [29L |30|. As it is evident, 
the two models can be easily differentiated from the dif- 
ferent g(z) behaviour, while their background expansion 
is almost the same as shown by the very close agreement 
between the corresponding p(z) and H(z) curves plotted 
in Fig. [3J According to the g(z) values, structures growth 
is faster in the tanh than in the exp model so that an ob- 
servational determination of g(z) can easily discriminate 
among the two f(T) scenarios. Indeed, using the data 
plotted, we get x 2 = 15.3 for the tanh model indicating 
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a strong disagreement, while better results are obtained 
for the exp model leading to \ 2 = 2.1 (still large, but 
mainly due to the highest z point which is not even fit- 
ted by the ACDM growth factor). These values have, 
however, to be taken with great caution and should not 
be considered as a definitive evidence of the exp model 
working better than the tanh one. First, most of the 
measurements of g(z) reported in [29j actually refers to 
g(z)/b(z) with b(z) the bias of the galaxy population used 
to trace the growth factor. Since the bias is related to the 
underlying gravity theory, a modified Lagrangian can in 
principle introduce deviations from the standard collapse 
scenario thus leading to a different and possibly scale de- 
pendent bias. Second, g{z) is sometimes obtained from 
a preliminary modelling of the matter power spectrum 
as in [30| which partially relies on the assumption of a 
fiducial ACDM model to convert from the redshift to the 
real space clustering. While the background expansion 
is likely the same as the ACDM one up to intermediate 
z, the modeling of the power spectrum and its distor- 
tions can be different so that the inferred g{z) estimate 
should be taken cum grano salis. Motivated by these 
considerations, we prefer to be conservative and consider 
Fig.[5]as a strong evidence that accelerating f(T) models 
with a quite similar background expansion can be easily 
discriminated by their growth factor evolution. 



V. CONCLUSIONS 

Interpreting gravitational interactions in terms of the 
torsion rather than the scalar curvature leads to the 
equivalent teleparallel formulation of GR. Adding a fur- 
ther f(T) term into the gravity Lagrangian then modifies 
General Relativity allowing for a wider range of possible 
dynamical behaviour. In particular, for some particular 
choices of the f(T) functional expression, it is possible to 
get an accelerating expansion in a matter only universe 
with the effective torsion fluid playing the role of dark 
energy. Motivated by these considerations, we have here 
considered two generalisations of teleparallel models as- 
signed by the expression (JBJ for the f(T) function. We 
have convincingly shown that both models are in excel- 
lent agreement with a wide set of cosmological data, from 
the SNela + GRB Hubble diagram to the H(z) measure- 
ments, BAOs and CMBR distance priors. 

Having been designed to give an accelerated expansion, 
the parameter space of both models collapses into a re- 
gion giving rise to a background dynamics similar to the 
ACDM one. As a consequence, the two models can be 
hardly discriminated based on the dataset we have used. 
It is worth stressing that this is not a limitation of the 
data, but rather an intrinsic feature of how the models 
have been worked out. As such, improving the precision 
of the measurements or increasing the statistics does not 
help in discriminating among the f(T) models and the 
ACDM one. On the contrary, one has to resort to dif- 
ferent tracers which are related to the evolution of the 



perturbations, the simplest one being the growth factor 
g{z). On small scales, the impact of torsion only intro- 
duces a redshift dependent rescaling of the gravitational 
constant which now becomes G e // = G/(l + fx)- This 
only modification considerably changes the growth factor 
leading to a scaling of g with z which strongly depend on 
the f(T) model considered. As Fig. [5] clearly shows, it 
is easy to discriminate between the two models provided 
one has a model independent estimate of the growth fac- 
tor as function of z. However, to this end, one should 
carefully take into account how the collapse of structures 
take place in the modified regime introduced by the ef- 
fective torsion field. To this end, the weak field limit of 
f(T) models must be investigated thus also leading to the 
determination of the gravitational potential on galactic 
scales and its impact on the galaxy dynamics. 



On the largest scales, the torsion field modifies the 
growth of perturbations by altering both the gravita- 
tional constant and the friction term [3(j. As a conse- 
quence a richer phenomenology is achieved possibly lead- 
ing to other ways to discriminate among f(T) models and 
dark energy ones. Two candidate probes are the mat- 
ter power spectrum P(k) and cosmic shear. In the first 
case, the derivation of P(k) is likely quite easy to achieve 
the main deviations from the standard GR one being en- 
coded by the scale dependent growth factor. However, 
the comparison with the data asks for two further effects 
to be taken into account. First, what we actually ob- 
serve is the galaxy power spectrum which is related to 
the matter one through the bias b (if we neglect redshift 
distortions). While for GR models one can assume b to 
be scale independent, this is not guaranteed a priori for 
f(T) models. As a possible way out, one could rely on 
a parameterized phenomenological ansatz although this 
strategy could introduce a degeneracy with the f(T) pa- 
rameters. On smaller scales, where g(z) is scale indepen- 
dent, nonlinear effects start becoming important so that 
one should find a way to take them into account. As far 
as we know it, a mapping from the linear P(k) to the non- 
linear one is still unavailable for f(T) theories thus pre- 
venting a straightforward comparison between predicted 
and observed P(k). A somewhat better situation holds 
for the computation of the cosmic shear power spectrum. 
Indeed, it probes the full matter distribution along the 
line of sight and not the clustered component. As such, 
the bias plays no role at all so that there is no need to 
preliminary investigate how galaxy formation takes place 
in f(T) models. Moreover, the cosmic shear power spec- 
trum allows us to probe the largest scales where nonlinear 
effects may be neglected thus putting more emphasis on 
the modifications due to the deviations of f(T) theory 
from GR. We therefore consider cosmic shear as a more 
promising tool to discriminate among rival teleparallel 
scenarios and between f(T) and GR and deserve a de- 
tailed analysis to a forthcoming publication |3lj . 
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